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摘 要 : 针对 陇 中 黄土 丘陵 沟 玖 区 土壤 水 蚀 过 程 复杂 且 难 以 有 效 预 测 的 问题 ,以 定西 市 安家 沟 水 土 


保持 试验 站 2005 一 2016 年 1 ~12 月 人 工 草 地 径流 场 试验 数据 为 主要 来 源 , 将 流域 月 降雨 量 、 月 侵 
蚀 性 降雨 量 、 月 径流 量 、 月 降雨 强度 、 径 流 场面 积 \ 径 流 场 坡度 、 土 壤 砂 粒 食量、 土壤 粘 粒 含量 8 个 区 


子 作为 输入 因子 ,月 土壤 水 蚀 量 作为 输出 ,运用 偏 最 小 二 乘法 (Partial Least-Squares Regression , 
PLSR) 和 长 短期 记忆 (Long Short-Term Memory ,LSTM) 循环 神经 网 络 建立 人 工 草 地 土壤 水 蚀 预 测 模 


型 ,并 利用 BP( Back Propagation ) .RNN( Recurrent Neural Network ) .LSTM 常见 神经 网 络 模型 ,对 模 
型 的 有 效 性 进行 评估 。 结 果 表 明 :PLSR 将 模型 8 个 输入 因子 减少 为 4 个 ,从 而 有 效 解决 LSTM 神 
经 网 络 模型 对 样本 数量 要 求 过 高 的 问题 ; PLSR 和 LSTM 神经 网 络 模型 的 结合 可 以 有 效 提高 模型 
对 人 工 草 地 土壤 水 蚀 过 程 的 预测 精度 和 收敛 速度 ,预测 结果 的 平均 相对 误差 小 于 49% ,相关 系数 高 
于 其 他 3 种 神经 网 络 模型 ,而 迭代 次 数 、 均 方 根 误差 和 平均 绝对 误差 均 低 于 其 他 3 种 模型 ;研究 发 
现 坡 度 对 人 工 草地 土壤 水 他 过 程 影响 较为 明显 ,降雨 量 小 于 25 mm 时 ,人 工 章 地 土壤 水 亿 量 不 会 
随 坡度 增加 而 明显 增长 ,但 当 降 雨量 超过 25 mm 时 ,人 工 草地 土壤 水 亿 量 会 随 坡度 明显 增加 。 
PLSR-LSTM 神经 网 络 土壤 水 蚀 预 测 模型 可 以 准确 预测 陇 中 黄土 丘陵 沟 敬 区 人 工 草 地 土壤 水 狂 量 ， 


为 该 地 区 水 土 流 失 的 准确 预报 提供 新 的 思路 和 方法 。 


关键 词 : 
文章 编号 : 


陇 中 黄土 丘陵 沟 确 区 是 我 国土 壤 水 蚀 最 严重 的 
地 区 之 一 ,土壤 水 蚀 造 成 该 地 区 土地 水 土 流失 严重 、 
沟 密 面积 增多 ,农田 遭受 破坏 ,严重 阻碍 了 当地 经 济 
的 可 持续 发 展 "。 研 究 表 明 降 雨 是 导致 坡 面 土壤 
水 刨 的 主导 因素 ,是 导致 水 土 流失 最 重要 的 外 部 驱 
动 , 而 坡度 、 坡 形 、 以 及 植被 覆盖 对 水 人 蚀 的 影响 也 极 
为 显著 ,整个 土壤 水 蚀 过 程 相对 复杂 ,难以 做 出 有 效 
的 预测 ”…… 。 目 前 ,国内 外 学 者 对 土壤 的 水 蚀 机 理 
和 规律 进行 了 大 量 的 研究 工作 ,并 在 此 基础 上 建立 
了 适用 于 不 同 地 区 的 土壤 水 蚀 模型 ,主要 包括 经 验 
统计 模型 和 物理 模型 ,其 中 经 验 统计 模型 多 以 多 元 
日 模型 为 主 ,模型 计算 简便 ,但 预测 精度 不 够 理 
想 ;而 物理 模型 的 建立 多 以 土壤 水 蚀 机 理 为 基础 , 模 
型 虽 能 较 好 的 反映 土壤 水 蚀 过 程 和 机 制 , 但 计算 较 
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近 几 年 , 随 着 深度 学 习 技术 的 不 断 成 熟 ,一 些 国 
内 外 研究 学 者 将 基于 深度 学 习 的 各 种 模型 应 用 于 时 
间 序 列 数 据 的 预测 中 "" 。 深 度 学 习 模型 是 一 种 由 
很 多 非 线性 隐藏 层 组 成 的 深度 神经 网 络 模型 ,能够 
根据 输入 数据 和 逐 层 提取 数据 的 有 效 特 征 , 从 而 挖掘 
出 数据 间 的 潜在 规律 。 循 环 神经 网 络 模型 ( Re- 
current Neural Network, RNN) 将 时 间 序 列 的 概念 引 
人 到 神经 网 络 模型 的 设计 中 ,使 其 更 适合 应 用 于 时 
间 序 列 数 据 的 处 理 "”| 。 在 深度 学 习 技术 不 断 发 展 
的 过 程 中 出 现 了 许多 RNN 模型 的 变 体 ,其 中 长 短期 
记忆 (Long Short-Term Memory, LSTM ) 神经 网 络 模 
型 通过 引入 “ 门 ” 解 决 了 RNN 长 期 记忆 能 力 不 足 的 
问题 ,使 LSTM 神经 网 络 模 型 真正 实现 了 长 距离 时 
间 序 列 数 据 的 有 效 利 用 。 但 模型 自 变量 间 的 多 
重 相 关 性 以 及 对 样本 点 个 数 要 求 过 高 等 问题 ,限制 
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区 人 工 草地 土壤 水 蚀 预 测 模型 


了 该 模型 的 实际 应 用 。 偏 最 小 二 乘法 ( Partial Least- 
Squares Regression, PLSR ) U“ 既 充分 吸收 主 成 分 分 
Br AR! ,同时 可 以 提高 主 成 分 之 间 的 相关 性 5 。 
利用 该 算法 可 以 实现 高 维 数据 空间 的 数据 降 维 , 特 
别 适合 在 数据 间 存 在 严重 多 重 相关 性 以 及 样本 点 个 
数 相对 较 少 的 模型 建 模 中 应 用 。 

本 文 则 在 构建 基于 PLSR 优化 后 的 LSTM 
(PLSR-LSTM ) 神 经 网 络 人 工 草 地 土壤 水 蚀 预 测 模 
型 ,以 实现 陇 中 黄土 丘陵 沟 餐 区 人 工 草 地 土壤 水 蚀 
的 准确 预报 ,为 该 地 区 水 土 流 失 的 准确 预测 提供 有 
效 的 技术 支撑 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

本 试验 将 径流 场 设置 在 甘肃 省 定西 市 安家 沟 流 
域 ,该 地 区 属于 典型 黄土 高 原 丘 陵 沟 蜜 地 貌 , 平 均 海 
拔高 度 2 250 m ,多 年 平均 降雨 量 为 385 mm ,年 均 蒸 
发 量 为 1 531 mm ,年 平均 气温 为 6.5 CC 。 该 地 区 土 
壤 类 型 以 灰 钙 土 为 主 ,土壤 母 质 为 黄土 ,质地 为 粉 砂 
SEA HE, AS , 极 易 发 生 土壤 水 刨 。 植 被 属于 
森林 草原 干草 原 区 , 人 工 草地 以 红豆 草 和 紫花 首 蒂 
1.2 试验 设计 

试验 所 用 径流 场 被 设置 在 安家 沟 流域 阴 坡 中 
部 ,面积 为 10 m x5 m, 坡 度 分 别 为 5°、10°、15°、 
20° ,每 个 坡度 径流 场 重复 2 次 。 径 流 场 0 ~ 10 cm 
土 层 土壤 的 砂粒 含量 和 粘 粒 含量 分 别 为 24.3% 和 
9.7% 。 人 工 草地 覆盖 植被 为 红豆 草 ,覆盖 度 达 到 
70% ~85% ,平均 为 79% 。 
1.3 试验 测定 项 目 和 分 析 方 法 

(1) 降雨 量 测定 

利用 流域 内 自动 气象 站 和 虹吸 式 雨 量 计 , 来 记 
录 每 次 降雨 的 雨量 .降雨 时 间 和 平均 雨 强 ,以 日 降雨 
HES 12 mm 的 标准 来 划分 水 蚀 性 降雨 量 ] 。 

(2) 径流 量 测定 

各 个 径流 场 通过 集 流 池 对 径流 产 沙 量 进行 观 
测 。 地 表 产 生 径 流 后 ,用 水 尺 读 取 集 流 池 内 的 水 次 ， 
再 乘 以 集 流 池 的 面积 ,计算 出 径流 总 量 。 

(3) 土壤 水 蚀 量 测 定 

将 集 流 池 中 泥水 搅拌 均匀 后 ,用 体积 为 500 mL 
的 取样 瓶 在 集 流 池 的 不 同 区 域 中 取水 样 3 次 ,对 取 
得 水 样 进行 过 滤 后 烘 干 称 重 , 测 得 水 样 的 泥 沙 含量 ， 


从 而 得 到 各 径流 场 的 产 沙 量 , 产 沙 量 与 面积 的 比值 
即 为 土壤 水 蚀 量 。 
1.4 偏 最 小 二 乘法 

偏 最 小 二 乘法 (PLSR ) 可 以 实现 多 对 多 线性 回 
归 建 模 ,适合 在 自 变 量 和 因 变 量 个 数 较 多 ,数据 间 存 
在 多 重 相关 性 ,以 及 样本 数量 较 少 的 情况 下 使 用 。 
一 般 情况 下 , 偏 最 小 二 乘法 与 主 成 分 分 析 法 一 样 在 
建 模 过 程 中 不 需要 选用 存在 的 个 成 分 4,ty，…， 
t,, ,通过 交叉 有 效 性 检验 提取 其 中 的 m 个 主要 成 分 
(<n) ,从 而 实现 性 能 较 好 预测 模型 的 建立 。 

交叉 有 效 性 检验 通过 每 次 舍 去 第 i 个 观测 值 
(i=1,2,3,…,n) FHA FAY n-1 个 观测 值 按 最 小 
二 乘法 进行 建 模 ,抽取 m 个 成 分 后 建立 回归 方程 ， 
然后 把 舍 去 的 第 i 个 观测 点 带 入 所 建立 的 回归 方 
程 ,得 到 y(j=1,2,…,p) 在 第 i 个 观测 点 上 的 观测 
值 j(m) ,对 i=1,2,3,…,n 重复 以 上 的 验证 , 即 得 
到 抽取 m 个 成 分 时 第 j 个 因 变 量 y (7 = 1,2, ++, p) 
的 残 差 平 方 和 "| 。 
PRESS,(m) = È lr; -为 0] = 1,2,…,p) (1) 


i= 


Y= (略为 为) 的 残 差 平 方 和 为 
PRESS(m) = YPRESS,(m) (2) 


当 PRESS(m) 达 到 最 小 时 ,对 应 的 m 即 为 所 求 
的 主 成 分 变量 个 数 , 而 本 研究 所 涉及 的 输出 变量 只 
有 一 个 ,所 以 令 上 式 p =1, 即 为 单 因 变 量 PLSR。 
1.5 LSTM 神经 网 络 模型 
LSTM 神经 网 络 模型 是 一 种 特殊 的 RNN 神经 
网 络 模型 ,通过 对 长 距离 时 间 序 列 数据 的 学 习 ,提高 
了 模型 预测 的 精度 。LSTM 神经 网 络 模型 通过 sig- 
moid 函数 和 一 个 点 乘 运 算 组 成 一 个 具有 信息 选择 
功能 的 门 ,其 中 sigmoid 函数 的 输出 值 在 0~1 之 间 ， 
0 代表 过 滤 信 息 ,1 代表 保留 信息 ,LSTM 神经 网 络 
模型 通过 门 能 实现 了 信息 的 遗忘 和 记忆 。LSTM 神 
经 网 络 模型 由 很 多 细胞 单元 组 成 ,每 个 细胞 单元 包 
括 遗 忘 门 .输入 门 和 输出 门 , 这 些 细胞 单元 使 LSTM 
神经 网 络 模型 具有 长 期 记忆 信息 的 能 力 呈 ,LSTM 
神经 网 络 模 型 细胞 单元 结构 如 图 1 所 示 , 其 计算 方 
法 可 表示 为 : 
i, =o (w, X[h,_,,x,] +b;) (3) 
f, =a (w, x [h,_,,,] +b) (4) 


tó wy 


c, =tanh(w, x [h,_,,x,] +b.) (5) 
c, =f, xc ti, Xe, (6) 
0, =a(w, x[h,_,,*,] +b,) (7) 
ae Sanh) (8) 


式 中 :x%, 为 本 细胞 单元 的 输入 值 ; 凡 -为 上 一 细胞 单 
元 的 输出 值 ;i, 为 输入 门 ;f. 为 遗忘 门 ;o, 为 输出 门 ; 
c, 为 LSTM 神经 网 络 模型 细胞 结构 状态 ;w 为 权重 
系数 和 矩阵;b 为 偏 置 项 ;o 代表 sigmoid PŽ; tanh 代 
表 双 曲 正 切 激活 函数 '” 。 该 算法 的 具体 计算 步 又 
为 : 

(1) 根据 上 一 细胞 单元 的 细胞 状态 c,_1 ,输出 
值 h_ ,以 及 本 细胞 的 输入 值 % ,按照 公式 (3) ~ 
(8) 计 算出 本 细胞 单元 的 输出 。 

(2) 按时 间 反 向 传播 法 ,利用 损失 函数 计算 每 
个 细胞 单元 的 反 向 误差 大 小 。 

(3) 根据 反 向 误差 ,计算 出 对 应 的 权重 梯度 。 

(4) 利用 基于 梯度 的 优化 算法 更 新 权重 降低 损 
PABA ,直到 LSTM 神经 网 络 模型 收敛 。 


h, 
Cm > C, 
f a 
o tanh 
h, I > 4 | > h, 
X, 
图 1 LSTM 神经 网 络 模型 细胞 单元 结构 图 


Fig.1 Cell structure of LSTM neural network model 


1.6 PLSR-LSTM 神经 网 络 模型 

首先 根据 PLSR 交叉 有 效 性 原则 从 输入 的 上 个 
自 变量 xz PRER m SEB ty yates 
经 “ 归 一 化 ”处理 后 输入 到 LSTM 神经 网 络 模型 的 隐 
WZ , RE ty ty yt ,i 输入 结束 后 ,LSTM 神经 网 
络 模型 的 细胞 单元 对 输入 数据 进行 计算 ,并 将 隐藏 
层 最 终 状 态 h, 传输 到 输出 层 , 经 “去 归 一 ”化 处 理 
后 得 到 最 终 的 预测 结果 y, 如 图 2 所 示 。 
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图 2 PLSR-LSTM 神经 网 络 模型 结构 
Fig.2 Structure of PLSR-LSTM neural network model 


2 结果 与 分 析 


2.1 PLSR 主 成 分 分 析 

本 文 将 径流 场所 在 流域 的 月 降雨 量 x 月 水 蚀 
性 降雨 量 、 月 降雨 强度 x; 月 径流 量 x eH 
BA xs \ 径 流 场 坡度 x6 \ 径 流 场 土壤 人 砂粒 含量 x 径流 
场 土 壤 粘 粒 含 量 xs 共 8 个 因子 作为 自 变量 ,月 水 蚀 
E y 作为 因 变 量 , 并 对 所 用 数据 进行 归 一 化 处 理 。 
根据 PLSR 的 交叉 有 效 性 原则 对 输入 变量 进行 主 成 


分 分 析 以 此 降低 由 于 影响 因子 共存 而 产生 的 相互 重 
复 信息 ,减少 输入 自 变 量 个 数 并 滤 去 噪音 。 

通过 径流 场 输入 变量 PLSR 主 成 分 分 析 可 以 看 
出 当主 成 分 变量 个 数 为 4 时 ,各 径流 场 的 残 差 平方 
和 PRESS 值 达到 最 小 ,如 图 3 所 示 。 根 据 交 又 有 效 
性 原则 ,本 研究 从 选取 的 8 个 自 变 量 中 提取 4 个 作 
为 主 成 分 输入 LSTM 神经 网 络 模型 。 

本 研究 以 月 水 蚀 性 降雨 量 x 月 降雨 强度 xy 
月 径流 量 x 、 径 流 场 坡度 x, NAZE, H kE y 
为 因 变 量 ,建立 基于 PLSR-LSTM 神 经 网 络 的 人 工 草 
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LL] (c) 径流 场 3 


残 差 平方 和 


I 2 3 4 5 6 F 8 
主 成 分 变量 个 数 
图 3 4 个 径流 场 不 同 成 分 对 应 的 PRESS 值 


Fig.3 Variation of PRESS against number of components in four runoff fields 


地 土壤 水 蚀 预 测 模型 。 


2.2 ”PLSR-LSTM 神经 网 络 模 型 的 训练 及 参数 
调整 


基于 PLSR-LSTM 神经 网 络 的 土壤 水 蚀 预 测 模 
型 将 定西 市 安家 沟 2005 一 2016 年 1~12 月 径流 场 
试验 数据 按照 训练 样本 占 70% ,测试 样本 占 30% 进 
行 划 分 ,其 中 2005 一 2013 年 1 ~ 12 月 试验 数据 作为 
训练 样本 对 PLSR-LSTM 神经 网 络 模 型 进行 训练 ， 
2014—2016 年 1 ~12 月 试验 数据 作为 测试 样本 用 
于 PLSR-LSTM 神经 网 络 模型 的 有 效 性 检验 。 

利用 训练 样本 对 PLSR-LSTM 神经 网 络 模型 进 
行 训 练 ,采用 正 交 化 方法 对 权重 系数 和 抢 阵 进行 初始 
化 ,将 月 水 蚀 性 降雨 量 x 、 月 降雨 强度 x 月 径流 量 
x, EMAI RE x, 输入 PLSR-LSTM 神经 网 络 模型 
得 到 月 土壤 水 蚀 量 的 预测 值 y' ,并 将 该 值 与 月 土壤 
水 蚀 量 实测 值 y 的 均 方 误差 作为 PLSR-LSTM 神经 
网 络 模型 的 损失 函数 ,损失 函数 如 下 所 示 : 


LL 起 ss 
loss = We (yi — Yi) (9) 
i=l 


DIRA PKZ loss 最 小 为 优化 目标 ,在 设 定 网 络 
初始 学 习 率 mg、 批 尺寸 以 及 最 大 迭代 次 数 step 的 条 
件 下 ,利用 Adam 下 降 优化 算法 经 反复 迭代 更 新 网 
RAN HE FA CAE BE , 降低 损失 函数 loss 的 值 , 直到 
PLSR-LSTM 神经 网 络 模型 收敛 ,PLSR-LSTM 神经 网 
络 模型 隐藏 层 参数 如 表 1 所 示 。 


ll [ (pb) 径 济 场 2 
10 | (b) 径流 场 
展 
B 
Be 
W 
x 
1 2; 3 4 5 6 7 8 
主 成 分 变量 个 数 
ll K @ 径流 场 4 
10 上 
9 bs 
8 L 
长 7| 
he 6p 
au 51| 
R 4 上 
3 L 
2 | 
1 1 1 1 1 1 1 J 
1 2 3 4 5 6 7 8 


主 成 分 变量 个 数 


表 1 隐藏 层 网 络 参数 


Tab.1 Network parameters of hidden layer 


参数 名 称 参数 值 
隐藏 层 数 3 

最 大 迭代 次 数 600 

批 尺寸 72 

损失 函数 均 方 误差 
下 降 方 式 Adam 优化 算法 
学 习 速 率 0.001 

权 值 初始 化 方式 正 交 化 初始 化 


2.3 PLSR-LSTM 神经 网 络 模型 的 有 效 性 与 复杂 
度 分 析 

为 了 验证 PLSR-LSTM 神经 网 络 模型 对 陇 中 黄 
土 据 陵 沟 蜜 区 人 工 草 地 土壤 水 刨 量 预测 的 有 效 性 ， 
本 研究 采用 均 方 根 误差 (RMSE ) 决定 系数 (及 ) 和 
平均 绝对 误差 (M4E) 3 种 统计 指标 来 进行 分 析 , 其 
中 RMSE „MAE 越 小 ,R 越 接近 1, 预 测 值 与 实测 值 
间 的 误差 越 小 ,模型 的 拟 合 效果 越 好 。 同 时 为 了 
分 析 PLSR-LSTM 神经 网 络 模 型 在 人 工 草地 土壤 水 
蚀 量 预 测 方面 所 具有 的 优势 ,本 研究 选用 BPO 、 
RNN 以 及 LSTM 等 常见 神经 网 络 模型 与 该 模型 的 
预测 结果 进行 对 比 ,结果 如 表 2 所 示 。 

PLSR-LSTM 神经 网 络 模型 对 4 个 径流 场 土壤 
水 蚀 量 预测 结果 的 RMSE 分 别 为 0. 234、0. 294、 
0. 321 .0. 368; MAE 分 别 为 0. 246 .0. 318、0. 348、 
0.381; R° 分 别 为 0. 985 .0. 982 0. 977 .0. 965, 3 种 
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R2 几 种 常见 神经 网 络 预测 模型 有 效 性 对 比 
Tab.2 Comparison of effectiveness of several common 


neural network prediction models 


径流 场 模型 RMSE MAE R? 
径流 场 1 PLSR-LSTM 0.234 0.246 0.985 
LSTM 0.354 0.371 0.973 
RNN 0.397 0.411 0. 958 
BP 0.424 0.457 0.943 
径流 场 2 PLSR-LSTM 0.294 0.318 0. 982 
LSTM 0.389 0.392 0.961 
RNN 0.418 0.439 0.944 
BP 0.465 0.493 0.935 
径流 场 3  PLSR-LSTM 0.321 0.348 0. 977 
LSTM 0.416 0.435 0.956 
RNN 0.442 0.481 0.941 
BP 0.493 0.524 0. 929 
径流 场 4  PLSR-LSTM 0.368 0.381 0.965 
LSTM 0.447 0.487 0. 940 
RNN 0.562 0.590 0.921 
BP 0.586 0.609 0. 906 


指标 值 均 优 于 其 他 3 种 神经 网 络 模型 的 预测 结果 。 
研究 结果 表明 ,PLSR-LSTM 神经 网 络 模型 能 够 有 效 
预测 陇 中 黄土 丘陵 沟 殉 区 人 工 草地 的 土壤 水 蚀 量 。 

为 了 对 PLSR-LSTM 神经 网 络 人 工 草地 土壤 水 
蚀 预 测 模型 的 复杂 度 进行 检验 ,本 研究 将 4 种 神经 
网 络 模型 的 损失 误差 进行 了 对 比 ,如 图 4 所 示 。 研 
究 结果 表明 ,4 种 模型 的 损失 误差 变化 趋势 相似 ,但 
PLSR-LSTM 神经 网 络 模型 的 损失 误差 小 于 其 他 3 
种 模型 , 量 该 模型 的 迭代 次 数 达 到 54 次 时 模型 的 损 
失误 差 趋 于 稳定 ,收敛 速度 高 于 其 他 3 种 神经 网 络 
模型 。 
2.4 PLSR-LSTM 神经 网 络 模型 土壤 水 蚀 量 预测 
分 析 

PLSR-LSTM 神经 网 络 模型 利用 2014 一 2016 年 
1~12 月 的 月 水 蚀 性 降雨 量 .月 降雨 强度 .月 径流 
Ht 径流 场 坡度 ,预测 当月 的 月 土壤 水 蚀 量 ,2014 一 
2016 年 月 平均 土壤 水 蚀 量 如 图 5 所 示 。 从 月 均 土 
坏 水 蚀 量 来 看 , 陇 中 黄土 丘陵 沟 完 区 人 工 草地 土壤 水 
刨 主要 集中 在 5 ~9 月 ,其 中 8 月 的 水 蚀 量 最 大 ,其 次 
是 7 月 和 9 月 ,4 个 坡度 径流 场 的 月 水 蚀 量 大 小 依次 
为 :径流 场 4> 径 流 场 3 > 径流 场 2 > 径流 场 1。 

将 月 土壤 水 蚀 量 求 和 转换 为 年 土壤 水 蚀 量 , 并 
计算 年 土壤 水 蚀 量 预测 值 与 实测 值 的 相对 误差 ,如 
表 3 所 示 。PLSR-LSTM 神经 网 络 模 型 预测 的 年 土 


壤 水 蚀 量 接近 于 实测 值 ,平均 相对 误差 小 于 4% 。 
由 此 可 见 , 将 PLSR 与 LSTM 神经 网 络 模 型 结合 可 以 
实现 陇 中 黄土 丘陵 沟 移 区 人 工 草地 土壤 水 蚀 量 的 准 
确 预 测 。 为 了 更 加 直观 显示 模型 的 预测 结果 ,将 
2014 一 2016 年 土壤 水 蚀 量 预测 结果 与 真实 值 进 行 
对 比如 图 6 所 示 。 从 年 土壤 水 蚀 量 大 小 来 看 ,径流 
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图 5 2014—2016 年 不 同 径流 场 月 均 土 壤 水 蚀 量 


Fig.5 Monthly average soil water erosion in different 


runoff fields from 2014 to 2016 
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Fig.6 Comparion chart of soil water erosion measured 


and predicted values from 2014 to 2016 
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#3 2014—2016 年 土壤 水 蚀 量 的 实测 值 与 预测 值 对 比 


Tab.3 Comparison of soil erosion measured and predicted values from 2014 to 2016 


径流 声 ani w 2016 年 

JAKE Hk AARE SAKE WKE ”相对 误差 AKE WAKE HIRA 
径流 场 1 67.8 66.4 2.06 18.6 18.8 1.08 26.2 25.8 1.53 
径流 场 2 100.8 98.6 2.18 27.6 27.2 1.45 38.8 37.6 3.09 
径流 场 3 176.8 182.8 3.39 56.2 55.2 1.78 112.6 109.2 3.02 
径流 场 4 215.8 225.6 4.54 167.0 172.2 3.13 238.4 249.8 4.78 
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图 7 不 同 降雨 量 条 件 下 人 工 草地 水 蚀 量 与 坡度 间 的 关系 


Fig.7 Relationship between water erosion and slope 


in artificial grassland with different rainfall 


场 4> 径 流 场 3 > 径流 场 2 > 径流 场 1 ,研究 结果 表 
明 ,坡度 对 土壤 水 蚀 的 影响 较为 明显 。 

不 同 降雨 量 条 件 下 土壤 水 蚀 量 与 坡度 的 关系 ， 
如 图 7 所 示 。 不 同 坡度 条 件 下 ,降雨 量 <10 mm 时 
引起 的 土壤 水 蚀 量 最 小 , 且 坡 度 增 加 对 土壤 水 蚀 量 
的 影响 不 大 ; “4 10 mm < 降雨 量 <25 mm 时 ,5° ~ 
20? 坡 度 区 间 的 土壤 水 蚀 量 随 着 坡度 的 增加 呈 小 幅 
度 增 加 ,增幅 并 不 明显 ; 当 25 mm < 降雨 量 <50 mm 
时 , 随 着 坡度 的 增加 ,土壤 水 蚀 量 增幅 加 大 ; 当 降 
雨量 >50 mm 时 ,土壤 水 蚀 量 随 坡 度 增 加 大 幅度 
增长 ,其 中 坡度 大 于 10°? 时 ,土壤 水 刨 量 的 增幅 最 
为 明显 。 


3 结论 


土壤 水 蚀 过 程 是 一 种 复杂 的 非 线性 过 程 , 由 于 
LSTM 神经 网 络 模型 具有 很 好 的 非 线性 映射 能 
因此 该 模型 可 以 被 应 用 于 土壤 水 蚀 量 预测 ,但 由 于 
模型 精度 对 学 习 速 率 取 值 敏 感 ,使 得 学 习 速 率 如 果 
设置 的 不 合理 可 能 造成 模型 预测 结果 出 现 从 拟 合 或 
过 拟 合 现象 ,影响 模型 对 土壤 水 蚀 预 测 的 准确 度 。 
将 PLSR 主 成 分 提取 与 LSTM 神经 网 络 模型 结合 ， 
分 发 挥 各 自 优点 ,通过 主 成 分 分 析 剔 除 观 测 数据 的 
重复 信息 ,在 有 效 降 低 数 据 间 相关 性 的 同时 ,也 解决 


了 LSTM 神经 网 络 模 型 对 实测 数据 要 求 过 高 的 问 
题 , 简 化 了 LSTM 神经 网 络 模型 的 系统 结构 ,提高 了 
模型 的 收敛 速度 和 可 靠 性 。PLSR-LSTM 神经 网 络 
模型 土壤 水 蚀 量 预测 结果 的 平均 相对 误差 小 于 
4% ,RMSE „MAE 和 R 等 统计 指标 均 明 显 优 于 BP. 
RNN LSTM 等 常见 神经 网 络 模型 。 研 究 结果 表明 ， 
该 模型 能 较 好 的 实现 陇 中 黄土 丘陵 沟 移 区 人 工 草 地 
土壤 水 蚀 的 准确 预测 。 

降雨 是 导致 陀 中 黄土 丘陵 沟 蜜 区 人 工 草 地 坡 面 
产生 土壤 水 蚀 的 主要 因素 之 一 ,同时 坡度 对 人 工 草 
地 土壤 水 蚀 的 影响 也 比较 明显 ,坡度 的 增加 会 导致 
土壤 水 蚀 量 的 增加 。 研 究 结果 表明 , 当 降 雨量 小 于 
25 mm 时 ,人 工 草 地 土壤 水 蚀 量 并 不 会 随 坡 度 增 加 
而 明显 增长 ,但 当 降 十 超过 25 mm 时 , 人 工 草地 土 
壤 水 蚀 量 会 随 坡度 明显 增加 。 

本 研究 所 采用 的 输入 变量 均 为 容易 获取 或 测定 
的 数据 ,但 在 预测 过 程 中 没有 考虑 植被 地形 等 其 他 
因素 对 土壤 水 蚀 预 测 的 影响 。 因 此 ,本 研究 采用 的 
数据 并 不 能 准确 地 还 原 土壤 水 蚀 形成 区 的 要 素 状 
况 ,使 得 土壤 水 蚀 预 测 过 程 存 在 误差 ,在 后 期 的 研究 
过 程 中 ,可 以 通过 增加 新 的 约束 条 件 来 进一步 提高 
模型 的 预测 精度 ,为 陀 中 黄土 丘陵 沟 蜜 区 土壤 水 蚀 
的 准确 预报 提供 有 效 的 技术 支持 。 


参考 文献 ( References ) 


[1] 王权 威 , 唐 莉 . 基于 ABC-BP 的 土壤 侵蚀 量 预 报 模型 研究 [J]. 
水 利 发 电 ,2017,43 (9) :1 - 4. [WANG Quanwei, TANG Li. 


Study on soil erosion prediction model based on artificial hee colo- 
ny algorithm and BP neural network [J]. Water Power, 2017, 43 
(9) 31-4. ] 

[2] 卫 伟 , 陈 利 顶 , 传 伯 杰 ,等 . 半 干 时 黄土 丘陵 沟 蜜 区 降水 特征 值 
和 下 垫 面 因 子 影响 下 的 水 土 流失 规律 [可 . 生态 学 报 ,2006 ,26 
(11) :3847 -3853. [ WEI Wei, CHEN Liding, FU Bojie, et al. 


Mechanism of soil and water loss under rainfall and earth surface 
characteristics in a semiarid loess hilly area [J]. Acta Ecologica 


Sinica ,2006 ,26 (11) :3847 - 3853. ] 


202006.00019v1 


chinaXiv 


[4] 


[6] 


[7] 


[8] 


[9] 


[13] 


tó wy 


张 兴 奇 , 顾 礼 彬 , 张 科 利 ,等 . 坡度 对 黔 西北 地 区 坡 面 产 流产 沙 
的 影响 [中 . 水土 保持 学 报 ,2015 ,29(4) :18 -22 ,72. [ ZHANG 
Xingqi,GU Libin,ZHANG Keli,et al. Impacts of slope gradient on 


runoff and sediment in northwest Guizhou[ J]. Journal of Soil and 
Water Conservation ,2015 ,29(4) :18 -22 ,72. ] 

和 继 军 , 蔡 强 国 , 刘 松 波 . 次 降雨 条 件 下 坡度 对 坡 面 产 流产 沙 
的 影响 [中 .应 用 生态 学 报 ,2012 ,23(5):1263 - 1268. [ HE Ji- 


jun, CAI Qiangguo ,LIU Songbo. Effects of slope gradient on slope 
runoff and sediment yield under different single rainfall conditions 
[J]. Chinese Journal of Applied Ecology ,2012,23 (5): 1263 - 
1268. ] 

ARAL, ME Bak, SE. BANAR rl ET ME 
[四 .水 土 保 持 学 报 ,2012 ,26 (4) :28 -36. [YANG Chunxia , XI- 
AO Peiqing,ZHEN Bin, et al. Characteristics on runoff and sedi- 


ment yield in field prototype and different coverage slope[ J ] . Jour- 
nal of Soil and Water Conservation ,2012 ,26(4) :28 - 36. | 

胡 琳 , 苏 静 , 桑 永 枝 , 等 . 陕西 省 降雨 侵蚀 力 时 空 分 布 特征 [了 . 
干旱 区 地 理 ,2014,37(6):1101 - 1107. [ HU Lin, SU Jing, 


SANG Yongzhi ,et al. Spatial and temporal characteristics of rain- 
fall erosivity in Shaanxi Province[ J]. Arid Land Geography ,2014, 
37(6); 1101 - 1107. ] 

EDL, Rie, BENE, E. TRMM 数据 在 中 国 降雨 侵蚀 力 计算 中 
的 应 用 [J 了 ]. 干旱 区 地 理 ,2015,38 (5 ) : 948 -959. [ WANG 
Kai, CHEN Lu,MA Jinhui, et al. Calculation of rainfall erosivity in 
China with TRMM data[ J]. Arid Land Geography ,2015 ,38 (5) : 
948 -959. ] 

RENARD K G,FOSTER G R,WEESIES G A,et al. RUSLE; Re- 


vised universal soil loss equation [ J]. Journal of Soil and Water 
Conservation ,1991 ,46(1) :30 -33. 

王 略 , 届 创 , 赵 国 栋 . 基于 中 国土 壤 流失 方程 模型 的 区 域 土 壤 
侵蚀 定量 评价 [J]. 水 土 保持 通报 ,2018,38 (1 ) :122 - 125. 
[WANG Lue, QU Chuang, ZHAO Guodong. Quantitative assess- 


ment of regional soil erosion based on Chinese soil loss equation 
model[ J]. Bulletin of Soil and Water Conservation ,2018 ,38 (1) : 
122 -125. ] 

LANGK V M,KARLSSON L, LOUTFI A, et al. A review of unsu- 
pervised feature learning and deep learning for time-series model- 
ing[ J]. Pattern Recognition Letters ,2014 ,42(1) :11 -24. 
LECUN Y,BENGIO Y, HINTON G. Deep learning [ J]. Nature, 
2015 ,521 (7553 ) :436 - 444. 

范 竣 翔 ,李琦 , 朱 亚 杰 ,等 . 基于 RNN 的 空气 污染 时 空 预报 模 
型 研究 [J] .测绘 科学 ,2017 ,42(7) :76 -83. [FAN Junxiang, LI 
Qi,ZHU Yajie ,et al. A spatio-temporal prediction framework for air 
pollution based on deep RNN[ J]. Science of Surveying and Map- 
ping 2017 ,42(7) :76 -83. ] 

EE, RIR, XH , S. ETF LSTM 循环 神经 网 络 的 故障 时 间 序 
列 预 测 [可 . 北京 航空 航天 大 学 学 报 ,2018 ,44 (4) :772 - 782. 


[14] 


[15 


i 


[16] 


[17] 


[18] 


[21] 


[22] 


ChinaXiv 合 作 期 刊 


[WANG Xin, WU Ji, LIU Chao, et al. Exploring LSTM based re- 
current neural network for failure time series prediction [ J ]. Jour- 
nal of Beijing University of Aeronautics and Astronautics 2018 ,44 
(4) :772 -782. | 

李 世 欣 , 温 建 , 邵 孝 侯 , 等 . 偏 最 小 二 乘法 与 人 工 神经 网 络 耦 合 
的 小 流域 产 沙 模型 [可 . 海河 大 学 学 报 ,2010 ,38(2) :149 - 
152, [LI Shixin, WEN Jian, SHAG Xiaohou et al. Sediment yield 


model for small watersheds based on coupling of partial least square 
regression and artificial neural network [ J]. Journal of Hohai Uni- 
versity ,2010 ,38 (2) :149 - 152. ] 

ei, NUTRI , XIVE , SE. 基于 PCA 和 BP 神经 网 络 的 径流 预测 
[J]. 中 国 沙漠 ,2016 ,36(4) :1144 - 1151. [NIE Min, LIU Zhi- 


hui, LIU Yang,et al. Runoff forecast based on principal component 
analysis and BP neural network[ J]. Journal of Desert Research , 
2016 ,36(4) :1144 - 1151. ] 
张 朋 , 下 崇 峰 , 杨 永 胜 ,等 . 基于 CCA 的 坡 面 尺度 生物 结 皮 空 
Hat LJ]. 生态 学 报 ,2015 ,35 ( 16 ) :5412 - 5420. | ZHANG 
Peng,BU Chongfeng,YANG Yongsheng ,et al. The slope scale dis- 


tribution regularity of biological soil crusts based on CCA[ J]. Acta 
Ecologica Sinica ,2015 ,35(16) :5412 - 5420. ] 

谢 云 , 刘 宝 元 , 章 文 波 . 侵蚀 性 降雨 标准 研究 [J]. 水 土 保持 学 
报 .2000,12(4) :6-11.[XIE Yun,LIU Baoyuan, ZHANG Wen- 


bo. Study on standard of erosive rainfall [ J]. Journal of Soil and 
Water Conservation ,2000 ,12(4) :6 -11. ] 

张 东 平 , 孙 克 勤 , 潘 效 军 , 等 . 烟 气 脱硫 系统 增 压 风机 仿真 模型 
发 及 应 用 [J]. 环境 工程 学 报 ,2011,5 (9 ) :2082 - 2086. 
[ ZHANG Dongping,SUN Keqin,PAN Xiaojun ,et al. Development 


and application of boost fan simulation model for flue gas desulfu- 
rization system[ J]. Chinese Journal of Environmental Engineering, 
2011 ,5(9) :2082 - 2086. ] 

GREFF K ,SRIVASTAVA R K, KOUTNIK J,et al. LSTM; A search 
space odyssey [J]. IEEE Transactions on Neural Networks & 
Learning Systems ,2016 , (99) :1 -11. 
王旭东 , 严 珂 , 陆 慧 娟 ,等 . LSTM 的 单 变量 短期 家 庭 电力 需求 
预测 [中 . 中 国 计 量 大 学 学 报 ,2018 ,29 (2) :142 - 147. [ WANG 
Xudong, YAN Ke, LU Huijuan, et al. Short-term household elec- 


tricity demand forecast based on LSTM single variable[ J]. Journal 
of China University of Metrology ,2018 ,29 (2) :142 - 147. ] 
ABUDU S,CUI C L,KING J P,et al. Comparison of performance of 
statistical models in forecasting monthly streamflow of Kizil River 
China[ J]. Water Science and Engineering , 2010, 3 (3) : 269 - 
281. 

刘 婷 , 邵 景 安 . 基于 BP 神经 网 络 的 三 峡 库 区 土壤 侵蚀 强度 模 
拟 [中 .自然 资源 学 报 ,2018,33(4): 669 - 683. [LIU Ting, 


SHAO Jing’an. Simulation of soil erosion intensity in the three gor- 
ges reservoir area using BP neural network [J]. Journal of Natural 


Resources ,2018 ,33(4) : 669 - 683. | 


F 


E PSE MPE ERR KATE ERK a RA 


Testing a soil water erosion predictive model in an artificial grassland 


in the hill-gully Loess Plateau 


WANG Jun'”, LI Guang', NIE Zhi-gang’, LIU Qiang’ 
(1 College of Forestry , Gansu Agricultural University , Lanzhou 730000 , Gansu , China; 
2 College of Information Science and Technology , Gansu Agricultural University , Lanzhou 730000 , Gansu , China ) 


Abstract; Erosion caused by soil water is widespread in the hill-gully Loess Plateau of central Gansu , China, but 
it is difficult to observe and predict. The objectives of this study were to establish a predictive model for soil water e- 
rosion based on a partial least squares algorithm (PLSR) and long short-term memory (LSTM ) , and to evaluate the 
model’s effectiveness using a back-propagation neural network ,a recurrent neural network ,and a LSTM neural net- 
work. The results were verified with experimental data from a runoff field of artificial grassland in the soil and water 
conservation experimental stations in Anjiagou, Dingxi City, from 2005 to 2016 (1 - 12 months). Input factors for 
the model were monthly precipitation , monthly erosive precipitation , monthly runoff, monthly rainfall intensity , runoff 
area ‚slope of runoff field, soil sand content, and clay sand content ,from which the model predicted monthly soil wa- 
ter erosion. The results revealed that the original eight input factors could be reduced to four by the PLSR algorithm 
in the prediction model, alleviating the requirement of the LSTM neural network model for a large number of sam- 
ples. Our method improved both the prediction precision and the convergence rate of the soil water erosion model in 
artificial grassland by combining a PLSR algorithm with a LSTM neural network. The average relative error of our 
prediction results was less than 4% , and the correlation coefficients were higher for our model than for the three 
neural network models used for comparison. In addition, iteration times , root mean square error, and mean absolute 
error were lower than for the other neural network models. We found that slope had a significant effect on soil water 
erosion , but only when the monthly rainfall was more than 25 mm. Our PLSR-LSTM neural network model could be 
used to accurately predict soil water erosion in artificial grassland in the hill-gully Loess Plateau of central Gansu, 
and it could provide a new path toward accurate prediction of soil erosion in this area. 

Key words; hill-gully loess plateau; artificial grassland; long short-term memory (LSTM) recurrent neural net- 


work; soil water erosion prediction 


